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This thesis is an examination of the short period fluctuations 
of tidal nature of the Florida Current surface flow from an analysis 
of direct surface current measurements made over a period of about 
a month in the Florida Straits during the latter portion of 1965. 
The principal conclusion reached is that the surface current is 
modulated by a diurnal standing wave, coupling the tides of the 
Atlantic Ocean with those of the Gulf of Mexico, which produces a 
pronounced diurnal effect on the surface current fluctuation. 
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lee INTRODUCTION 


The pulsation of marine currents has long been a topic of major 
interest in the field of physical oceanography. Much of the dynamics 
of oceanic circulation can only be understood in terms of the velocity 
field and the forces acting upon it. In general, the study of current 
fluctuations is best considered on the basis of information concerning 
current velocities. Unfortunately, reliable information of this type 
has not been available in the past because of the inherent difficulties 
and expense involved in the gathering of such data. Numerous analyses 
of oceanic current fluctuations have of necessity been made on data 
Over short duration periods, non-synoptic surveys, or resorted to the 
indirect measurement from dynamic computations. Conclusions based 
upon such information have often been highly speculative in regard to 
current fluctuations, although in many cases commendably accurate. 

Improved technology in the field of ocean measurement instrumenta- 
tion has made possible the collection of data over sufficient time in- 
tervals to provide meaningful results. The development of moored buoy 
systems in conjunction with dependable current recorders has provided 
a powerful tool which now permits us to commence directly examining 
the variation of ocean currents over long periods of time. The purpose 
of this paper is to examine the nature of current velocity fluctuations 
in the Florida Current using current speed measurements collected over 
a period of about a month at hourly intervals in the Straits of Florida 


from a moored buoy. Short period variations of tidal period will be 


investigated as a step towards the ultimate goal of long period fluc- 
tuation studies from which perhaps someday accurate predictions and 
forecasting of the Gulf Stream System and its associated effect on 
climate may be gained. It is evident that our knowledge of long- 
period variations will remain unsatisfactory until the short-period 


oscillations have been sufficiently described and explained. 





a le HISTORY 


The pulsations of tidal origin of the Florida Current have long 
been a subject of interest and speculation. As early as the mid 1800's, 
Pillsbury on BLAKE undertook the admirable task of occupying six an- 
chor stations between Fowey Rocks, Florida and Gun Cay, Bahama 
Islands, devoting over 1100 hours in observing the Florida Current at 
this section in the course of two summers (Pillsbury, 1891). The 
longest continuous anchorage was for a period of 166 hours. He found 
monthly variations related to the declination of the moon and daily 
oscillations of tidal origin in the Florida Current, the latter 
amounting in some instances to a variation as large as 128 cm/sec. 

He reported two periods of increase and two periods of decrease in the 
Florida Current speed during the period of a lunar day. 

In April 1937, Parr anchored ATLANTIS for a series of five hydro- 
graphic stations which were successively occupied for twenty-four hour 
periods in the Straits of Florida between Fowey Rocks and Gun Cay (Parr, 
1937). He reported from analyzing hourly surface current speeds that 
there was up to a 50 cm/sec fluctuation of the speed of the Florida 
Current that was apparently caused by tidal forces. He found both 
diurnal and semidiurnal variations with the latter predominant. 

From May 1950 to May 1951 the U.S. Navy Hydrographic Office 
recorded hourly Loran fixes from oil tankers traveling between Cape 
Hatteras and the Key West vicinity. From an analysis of over 5000 


observations, it was concluded that there was an apparent tidal perio- 





dicity in the surface current flow (O'Hare et al, 1953). An attempt 
was made to screen the tanker survey data for a semidiurnal tidal 
periodicity, but the results obtained were inconclusive. 

In 1951 Murray attempted to analyze velocity fluctuations in the 
Florida Current by making numerous crossings between Miami and Gun Cay 
utilizing the Geomagnetic Electrokinetograph (GEK) (Murray, 1952). He 
was unable to detect any apparent periodicity in the fluctuations of 
the Florida Current, but pointed out the inherent difficulties involved 
with attempting to observe periodicities from scattered GEK observations. 

From December 1952 to November 1953, GEK fix stations were oc- 
cupied between the Miami Sea Buoy and Gun Cay. An attempt was made to 
study tidal fluctuations in the Florida Current by plotting 18/7 GEK 
measurements against time (Hela and Wagner, 1953). The data seemed 
to indicate the existence of a tidal fluctuation, although the fluc- 
tuations were strongly masked by non-tidal effects. 

From studies of electric potential measurements between Key West, 
Florida and Havana, Cuba, Wertheim was able to show evidence of the 
diurnal tidal influence in the transport through the Florida Straits 
(Wertheim, 1954). He found that the ratio of the amplitudes of the 
Reammenic coefficients for the tidal components My, 59, K , and 0, were 
in between those of the semidiurnal Atlantic tide at Miami and the 
diurnal Gulf of Mexico tide at Galveston, confirming the dependence of 
the transport on both tidal systems. 

Webster analyzed a large amount of GEK data for both the Straits 
of Florida and off Onslow Bay, North Carolina. He concluded that 
although it was probably rash to ascribe the velocity fluctuation of 


the Florida Current predominantly to tidal causes, the periods of the 





fluctuations observed were in the order of one day (Webster, 1961). 
Recently Schmitz and Richardson utilized a least square harmonic 
analysis on transport data acquired over a period of three years using 
the free instrument technique across a section of the Florida Current 
(Schmitz and Richardson, 1967). Based on the limited data available, 
they indicate that it is possible that fluctuations of tidal period 
are the major modulation of the Florida Current transport. Their 
estimates of tidal coefficients for transport amplitudes are 3.5 + l 


(10°M.3 Sec”) for My, 0, and Kj, and 1.5 (10°M.3 sec”") for Sp. 





IIE EXPERIMENTAL PROCEDURE 


GENERAL 

During the latter portion of 1965 a General Dynamics ocean buoy 
was moored for testing and evaluation in the Straits of Florida at a 
location as shown in Figure 1. From the period November 20 to December 
18, 1965, the buoy was equipped with a rotary current meter placed inm- 
mediately beneath the water surface, installed by Dr. William 
Richardson. Current speed data for one minute averages taken ap- 
proximately each hour were telemetered to the mainland and recorded. 
Also available for this period are wind speed, wind direction, baro- 
metric pressure and the significant wave height at the time of re- 
cording. Current direction is considered as essentially steady in 
the North-South orientation of the Florida Current at this location. 
A plot of current velocity versus time is illustrated by Figure 2. 

Cursory examination of Figure 2 reveals fluctuations of current 
speed in the order of 30 cm/sec with apparent periodicity. Since the 
periodicity appears to be largely on the order of 24 hours and can be 
initially presumed to be due to tidal effects, it was decided to con- 
duct a harmonic analysis on the data to determine if the fluctuations 
were indeed due largely to the tidal effects. 
A. HARMONIC ANALYSIS OF TIDAL CURRENTS 

Harmonic analysis is a method for describing periodic phenomena 
in which values of the dependent variable repeat themselves at equal 


intervals of time. The first practical application of harmonic analysis 








Figure l. General Dynamics Monster Buoy Location 





















































= = — ER a = —— a eer = 
—— = = ; : W aa é Al } 
“44 i be ee 
J a , I io ee a Poo ee —_—_—o a “ws 7, eres ee 
_ —_— i es ee TT we ae eG i 7 al 20 | 
 ——— Es 7 = i ee ee Se 
-- Rd —_— = eae ee ee ee = 7 oe tT aj s u < 
'* Ff 436 ace a eee th vos a ) 
"te ul ee . Z 
a ‘ we 
> é ib Pan i = pile Ize . "GS e t y y 6 
HS Sases 249 acs 260 ? eee | 2” ay 
sd u Ps ¥ e t we. ‘ 
Ss 3 bis the ve a 
re nd \ Weeal of o Gt « 
= pil a re ‘ P 303 (Fa de \ $ wy ‘ | 
P C 33 a4 oe Maw, 
,** | : — Indien Pea wl ~) 
: 2 5 227) ee 
‘ wl! 4 , Betilement Pt’ “os 
a #36 324 FI 4 sec 44 fh GFA is 
: ae eee te 178 } aoe MAS 5 
4s 93 o 27 : FR & FIA 
=0% 
~ fn} 132 ae en 
aad { 302 | 
ys a f a3 : 
8 
as ews f aii 
# ahi: ‘get 213 rey) 
35) 
¢) — a 338 ' 
. 106 ; 435 ; 
a tsa | ats ‘ 
Sa? | 
— 
j fs ‘o *~ 235 D2 a ee 
' ies | ait ave 
-“ 4 320 209 
* 7402 : 
ae 233 | 269 | TT 345 
i 
| E 6 
t a &i2 | a° 264 
| Ler i t 
gai 
a H 200 
“3 24 a: 287 | 
a 4 4 20 
410 { _ 
¢2 13, | — \ m4 | ‘ 
4 276 | 333 c bi 
== q 
qa! i 
A 232 Wi set 
sul 313 | 225 350 | 2» os Fy b. if 
i | ri wae [ ~ ie 
“ | a ait : 1% & ia 
‘ 192 ae 183 “l| i 
’ t 
4 170 = 10 if 
27S a5 21% 19@ y : i i 
230° | Q 
! 251 by 
Ws Od a't = 
| ! : 
287 180 
rT: ’ 266 
145 J | 8 
i98 260 . as 
ae ra | « bing 
Looe ole ; ir lle +) ny ae 2, Q 
. | MONSLER & 10 al Hae : 
ie tre = 
o> 3s We Aw) Last 2 ) 210 . ge 
* 283 : ae Oe: CY ea ae 
es “— ut, <= 5 “A + 12 ’ t 5 ce aad « 
yo ‘ ’ Ly 460 <0 Bt ——- 4 5 
te3 ‘ \ L A i" fi) ( » : . [oe 2 ee a z 5 ne " 
> im cow ¥w. al ny ah & 
» ae » L. é “(> c » = me a : ae oe 4 \ eee 
2] 230 5 os eo 
x a 4 ai a aj 
Ze) . a \ ct 
cf 402 320 ** Saar ab 3k yee 
~ ay 3 422 a " { 
3} Séu ars “1 ‘ | a oo | 
» So; 346 3 at 3 mat 
3 : r) wa 
nz G4 Nae ey mI 
133 a 280 aa 3 . — . ai | 
3] 450 | x 7 aa 2! Hat 
e sow i" &G 4 Uae 
{ | 197,* piesa ae ae “ a . \. 4 , 
0 >| 490 427 fe Ga { } Wb eT \ \ rf | 
2° % ' ‘s Se Shan x : 
x Fl veee att AM eb ve en : 2 ‘ 7 | Cw sj, 
{ 0 = 290 200 / } 4 1 2 » +4 | 
304 - wo t | ae a a ra Ins 
s 479 : 2 ‘ rr Ct \ 21 oh? Ht 
103 5 North Darna 2 Bes 
109. ’ d Gs “ely, 
320 fs? id ee ee . . 
7 ram che > et od \ ' iT} ay + 34 
-" ; | wale eA lGn FUR (=+)20s06 281f8 23M | 
£ 490 in 408 it fee ROG Bn 396 Bai" a3 ae 
a i : Cana 
| 425 , at ' south Bimini "4 FI ' i / su 
35 i 410 BGS eee : i @ ' o th Sa 
18@ 420 : sos Mande as 2h ty | 
5 om A Pri yee | 2 | 
od a6 24 :| 232 t Mow ke i ' li Cis res | 
193 aA i S 410 440 4359 “ Aer WR ast 2 i td 4 
83 4 a H 440 BT bebe Cor | ' ti 2 "2 | 
| | i. = Set et Ty J 
1%) 106, Dae : ih &¢ 7 
eee ‘> t : = \e Il, ey VOrec eott lore i a ‘ mb 
L he y ee fu oT a. Bis ee et Cor ; . ' 
; : 4 ; 431 hale ae - | i‘ : 
on 2409 1> F 2.9 ES 443 Momm bers va "1 4 af ese is : y 
a ardind Kook 
\ 2 : 450 | er a 12 ui ua sin 
4 P ry | = \ vi or ' t . iF eee 7s bow 
7) / 2%) han iciee? Cayeji2; aa | } 
mn 











2 ee 2 ™ 
3 \ ' rT} 2 Re 2 I 
ne \ 2 2 atu ‘ fr 
Trumyh Rat / é‘ rey "480 9 | ay ia af af 
287 “he 2 \ fu Sg Fel 
ian ‘ \ 52 iy 34° ar Ces 12) 2 20 24 re | 
« ! 346 ; ie- 2 a 2 a 
y ; 9, A ition net ay ce ee ee oo ee ws 
457 © red) | 
RJiw 4 a. aN a — ‘7 
PL / 093 Ne oe, 0 © - ht a & 44, 
my Pnrets Hoel 177 Vf 449 pt ; aaa oat ER ans 20: t 
ee —— es See | ar aa 
ai 0 -rico To pe Te, Mita Sy ee my 
— a Cry “ ay’ cial 
a 50" 40 (CONTINUED 30" ON CHART buy ok 
w : 








a 
- —_—_- 
— 

a - en 





S eee oe 
ae + eee 
a 


a a 
eS Te et Se. 
- - ee 
LES — 
= 








Figure 2. Current Velocity versus Time for the Period 
21 November through 17 December 1965 in the 


Florida Current. 
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Figure 2. Current Velocity versus Time 





is attributed to Lord Kelvin who in 1867 devised a method of reducing 
tidal height observations to harmonic constituents (Darwin, 1898). 

The harmonic analysis of tidal currents or of currents having 
tidal components is a process by which the tidal components having 
relations to astronomical conditions are separated into elementary 
harmonic components or constituents. Each constituent represents a 
cyclical change during a particular period calculated from astronomi- 
cal data. The constituents are in reality a substitution of hypotheti- 
cal tide-producing satellites (having either fixed circular or ellip- 
tical orbits around the earth parallel to the equator) for the actual 
tide-producers, the moon and the sun. 

Theoretically, there are a large number of tidal constituents 
needed to accurately resolve the complicated motions of the moon and 
the sun into simple components. Generally, however, in any given 
location most of these are of small amplitude, and all but a few may 
be disregarded for practical purposes. The major constituents used 
to determine the principal features of tidal currents are listed in 
Table 1 along with their respective periods. Although the amplitudes 
and phases of the tidal constituents are modified by local causes, the 
period of the constituents is defined by the configuration and periodic 
relative motions of the earth-sun-moon system, and is invariable. In 
order to determine the tidal constituents listed in Table 1, observa- 
tions over a complete period of about fifteen days are required to 
Separate the individual diurnal and semidiurnal constituents with a 
fair degree of accuracy. It has been proposed that extremely accurate 
measurements coupled with an extremely low noise level would permit the 


analysis of shorter periods (Munk .and Hasselman, 1965), but such data 
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is not as yet available for the Florida Current. 

The harmonic constants sought for each of the major tidal con- 
stituents are amplitude and phase lag. Once determined, the constit- 
uents are then recombined and the resultant is the predicted tidal 
component of the current. 

Mathematically in tidal current analysis, a finite trigonometric 
sum is fitted to a set of current data in which there is a partial 
tidal periodicity. A finite number of points will exist since the 
current data will exist only at discrete points. Therefore, the data 
may be fitted by a finite number of sine and cosine curves. Solving 


for the five constituents in Table l this can be represented by 


> 
Lae) = ee + i [a,Sin Wyttb,Cosw,t | + Re (1) 


where 


Z(t) = periodic variate with time. (eg. current velocity), 


a ) = mean value of the variate, 

Oo 

n = the identification of the tidal constituent in equation (1), 
; th _ Fr 

W, = the angular velocity of the n constituent W, = = 


T 


where T is the period of the constituent), 


t = time, 
a and oe = Fourier coefficients for each constituent, 
and RR, = the residual 


The harmonic Fourier coefficients for the series are solved from 


9 N 
a =a Zu) Sin W ty (2) 
k=] 
2 N 
and lo. aaa APS COS Wy th (3) 
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where 


ct 
] 


k discrete points in time at which the data was collected, 


o 
ms 
Q. 
Zm 
i 


the total number of data points. 


It is convenient to combine the sines and cosines of equation (1) 


which belong to the same nth constituent (Gerges, 1966) into a single 
term 
A, Sin (Wh Ce I On) 
yielding 
5 
Zt) =a =S A, Sin (W,t, + @n) (4) 
n= 
where l 
_ 2 Jee, 
= [ae be | 
and On = == 


—arctan be ‘ 


A, is called the amplitude and p the phase of the ne constituent. 


In tidal analysis, A, of the major tidal constituents is cor- 
rected by a node factor which compensates for the effects of the 
variation in the moon's node in each 18.6 year cycle, and is also cor- 
rected for interference effects in the analysis from other constituents. 
These include significant lesser constituents which sometimes may not 
be resolved from the actual data because of the length of the data 
series, but which may be inferred from the solved constituents. This 
is permitted by the fact that although the amplitude and phases of the 
tidal constituents at any location differ considerably from their 
theoretical equilibrium theory values, amplitudes and phases of the 
constituents of nearby frequencies at any place have relations that, 
in general, agree fairly closely with the relations of their theoreti- 


cal coefficients. (Schureman, 1941). 
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 . is corrected for equilibrium arguments for starting times 
(to Zeta) in accordance with astronomical stages, and likewise cor- 
rected for interference effects from other tidal constituents to 
(Kappa). 

There are a variety of methods available for solving for the 
tidal constituents. In the classical form of tidal analysis, a series 
of current observations is divided into periods equal to the known 
period of each individual tidal constituent. Each period is then 
further sub-divided into a number of equal parts called the constituent 
hours which are numbered consecutively starting with zero at the 
initial instant for each period. Thus, the phase of each constituent 
and its overtides (harmonics) will be the same for a fixed instant in 
all sub-divisions having the same number, but all other constituents 
will have different phases (Dronkers, 1961). By summing and averaging 
all velocities which are observed at a fixed instant in each sub- 
interval (constituent hour) with the same number, the effect of all 
other components having incommensurable speeds as well as random noise 
will be predominantly eliminated. Convenient stencils are available 
to perform the separation of the various constituents from the current 
data recorded on proper forms. Each constituent is examined separately, 
and then the interference effects from other constituents are calculated 
and used to modify the results. A detailed description of the mechanics 
of this approach may be found in Schureman, 1941. This method of hand 
analysis becomes quite laborious and time consuming for all but the 
shortest of time series, but for over fifty years had been the estab- 
lished method of determinir tidal constituents. 


The classical approach has been computerized in recent years with 
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the only major improvement (other than in computing speed) being that 
successive equally spaced data is multiplied by sines and cosines of 
angles incremented by the exact angular speed of each constituent, 
whereas with the hand analysis, successive values are allocated by 
stencil to the nearest constituent hour and then later corrected for 
this approximation by the application of suitable augmenting factors. 

Another modern technique has now largely been adopted which is 
known as the least-squares method of tidal analysis. Such an approach 
utilizes a program which specifies the frequencies of the constituents 
which are being sought, and the harmonic coefficients of the tidal 
constituents are then determined simultaneously so that the average 
values of the residuals squared is a minimum. Correcting the results 
for node factors or the equilibrium arguments for the starting times 
remains the same as in the traditional approach. 

Comparative tests show that the harmonic constants for the same 
set of constituents are slightly more accurate utilizing the least- 
Squares method as opposed to the classical approach (Zetler and Lennon, 
1967). The greatest advantage in using the least-squares method, how- 
ever, is that it requires neither equally spaced data nor a synodic 
period, whereas the traditional Fourier tidal analysis requires both 
equally spaced data and a quasi-synodic period of the principal con- 
stituents (Zetler et al, 1965). In as much as in oceanographic data 
collection these latter conditions are usually difficult and sometimes 
impossible to fulfill, the least-squares method of analysis is often 
the only method which may be used to determine the tidal constituents. 
The least-squares method has now been adopted by the Coast and Geodetic 


Survey for the analysis of data series of one year's duration and is 
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sometimes used with shorter series. 
B. TOTAL CURRENT FLUCTUATION FIELD 


Over a Series of current measurements, the mean value of current 


speed is given by: v= KV/n 
“Vv = mean current speed 
v- = speed at any observation 
N = total number of observations 


The degree to which numerical data tend to spread about a mean 
value is termed the variation or dispersion of the data. One measure 
of the dispersion of the data is the standard deviation, s, which is 
often termed the root mean square deviation and which is defined by 

g¢ = Tir - 7)? 
N 

The variance, which is another method of measuring the spread or 

dispersion of the variate with reference to the mean, is defined as 


the square of the deviation or 5? 


» and the energy of all of the fluc- 
tuations of the current over a series of time is equal to the variance. 
The energy contributed by any periodic component of the fluctuation 


field is equal to 1/2 (a," + eo wiere a2. and b, are the Fourier 


n 
coefficients of any periodic component over the period of investigation. 
(For derivation, see Appendix A). 

Parseval's Theorem tells us that the energy of a composite wave 
is composed of the sum of the energies of each of the distinct harmonic 
constituents of the waves of different frequencies making up the basic 
wave (Ippen, 1966). Thus, from any current fluctuation field, the per- 
centage of total energy due to periodic fluctuations may be calculated 


and compared to the total energy of all of the fluctuations, and the 


total energy contributed by individual periodic components may be 
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IV. RESULTS 


A plot of current velocity versus time for the period 21 
November through 17 December 1965 obtained from the General Dynamics 
Monster Buoy in the Florida Straits was presented in Figure 2. Figure 
3 presents a plot of current velocity versus time with a superimposed 
plot of the predicted tidal component of the current from the results 
of this analysis. The predicted tidal current is plotted about the 
mean current (167.05 cm/sec) indicating the predicted Florida Current 
in the absence of fluctuations of a non-tidal nature. Figure 4 is a 
plot of the residual after the predicted tidal component of the fluc- 
tuations of the Florida Current has been removed. Hourly values for 
current speed, wind speed and wind direction for the 15 day period are 


tabulated in Appendix B. 


A. DISCUSSION 

As can be seen from Figure 2, the second half of the current 
data series is of much poorer quality than the first half. Gaps in 
the record were due to telemetry and recording failure. The intense 
abrupt fluctuations seem to be indicative of a gradual failure of the 
current meter. Data obtained from the first fifteen days of this period 
is continuous, and the current meter one minute speed averages were 
available for each hour. 

An initial computer analysis of the entire data series utilizing 


a least-squares program indicated that the diurnal tidal constituents, 











Figure 3. 


Current Velocity versus Time for the Period 21 
November through 5 December 1965 with Super- 


imposed Plot of the Predicted Tidal Modulation. 
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Figure 3. Current Velocity versus Time with Superimposed 
Plot of Predicted Tidal Modulation 
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Figure 4, 
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Plot of the Residual After Tidal Modulation is Removed. 
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Table 2. 


Comparison of Least-Squares Harmonic Analysis 
of the First Half of the Data Series with That 


of the Entire Series. 
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TABLE 7 


Comparison of least-squares harmonic analysis of the first 
half of the data series with that of the entire series 


1. 15 DAY ANALYSIS KL 


One 
(360 Data Entries) Amplitude* Phase Amplitude* Phase 
Og 59.8° O15 TUS ge fa 
2. ENTIRE DATA SERIES KL Ol 
(596 Data Entries) OS 2 a Oats 296.5° 


* Amplitude is in knots. 


Amplitudes and phases are values prior to correcting for 
equilibrium argument, interference effects, etc. These 
comparative results will slightly differ from those found 
elsewhere in this paper as during the above analysis the 
constituent Ny was included in the calculations. 
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Kj] and 0], are the predominant contributors to the tidal fluctuations 
of the current. The results of the least-squares analysis for the 
diurnal constituents for the first fifteen days and the entire period 
are compared in Table 2. That the phase angles change comparatively 
little between the two separate analyses is indicative that the 
amplitudes of the diurnal tidal components are in fact real and not 

due to random fluctuations. Because the second half of the data series 
was of marginal reliability, it is not included in the remainder of 

the analysis. 

Table 3 indicates that three techniques for solving for the 
major tidal constituents produce, as previously discussed, relatively 
small differences. These differences are accounted for in the method 
of each technique. The Fourier computer analysis was selected for 
determining the harmonic constants for the data series. The final 
results are contained in Table 4. A total of 24 constituents, of which 
20 were inferred from the major constituents, were then recombined to 
produce the predicted tidal modulation of the Florida Current for the 
period 21 November through 5 December 1965. The predicted plot of 
tidal modulation is shown superimposed on the basic current data in 
Figure 3. 

Astronomical data for the period is given in Table 5. It can be 
seen that the predicted tidal modulation is in agreement with the astro- 
nomical data for the period. The additive effect of My and S, on pap 
November is not readily apparent due to the semi-diurnal components 
being over-shadowed by Kj, and 0, as they approach phase agreement on 
26 November. The predominant diurnal modulation becomes negligible as 


K, and 01 come into opposition on 3 December leaving only semi-diurnal 











Tabwe 3. 


Comparison of Results of Harmonic Analysis on the 
First 15 Day Data Period Using Different Analysis 


Methods. 
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Table 4. 


Harmonic Constants from Florida Current Data. 
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TABLE 4 


Harmonic constants from Florida current data 


CONSTITUENT** AMPLITUDE (H) (cm/sec) KAPPA**«%* 


(a.) DIURNAL CONSTITUENTS 


Kl 5.710 24.49° 
01 5.551 12.25° 
Py 1.888 24.49° 
Qi 1.075 erie 
J] 0.437 30.61° 
M1 0.396 esse a 
001 0.237 36.74° 
RHO} 0.211 6.98° 
20; 0.144 0.00° 

(b.) OTHER CONSTITUENTS 
M2 3.400 oS 
$9 2.500 309.40° 
S¢ ih Wei 159.88° 
Ko 0.679 SO0e Ur 
Mg 0.664 277.76° 
No 0.658 203.68° 
Sh 0.458 236523" 
My, 0.304 ANEW Re 
Mg 0.278 290.55° 
T) 0.149 309. 40° 
NU> 0.129 eS ek 
Lo 0.098 26022 
2N> 0.087 292.41° 
LAMBDA? 0.026 Bomise 
Ro 0.020 309. 40° 





«xk Nomenclature and constituent speeds are in accordance with the 
classical Doodson classification. 


**x* To convert to Greenwich Epoch (G), add the product of 79.85 times 
the value of the constituent subscript. 
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Table 5. Astronomical Data for the Period 21 November 


Through 17 December 1965. 








TABLE 5 


Astronomical data for 
period 21 November 
through 17 December 1965* 


22 November New Moon 

26 November Moon farthest South of Equator 
29 November Moon in apogee 

1 December Moon in First Quarter 

3 December Moon at Equator 

8 December Full Moon 

10 December Moon farthest North of Equator 
11 December Moon in perigee 

15 December Moon in Last Quarter 

16 December Moon at Equator 


*« From Americen Ephemeris and Nautical Almanac. 
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tidal modulation of lesser than usual amplitude as Mj and S») were in 
phase opposition on Ll December. 

Fluctuations of the observed and predicted current are in excel- 
lent agreement and clearly show that the tidal modulation is the major 
periodic fluctuation of the current during the period 23 through 30 
November. During the remainder of the data series correlation with 
the predicted tidal modulation still remains fairly good, but the mean 
current speed is raised or lowered apparently in response to local 
wind stress or possibly in response to atmospheric variations over the 
water regions which are coupled to the Florida Straits. 

Figure 4 is a plot of the residual after the predicted tidal 
modulation has been removed from the data. Included is the mean local 
wind speed and average wind direction at the buoy location in the 
Straits during each day from midnight to midnight. It would appear 
that response of the current to wind stress from either the East or 
West is negligible, and further that the current fairly rapidly re- 
sponds with increased velocity to winds from southerly directions. 
Response of the current to northerly winds appears more complicated 
with the indication from this limited data series being that the cur- 
rent is slow to initially respond to northerly winds, but once the 
response commences the current speed is greatly lowered and recovery 
to normal conditions is quite gradual. The seemingly erratic period, 
30 November through 3 December, contained the highest wind speeds and 
rough seas were prevalent. 

Table 6 presents statistical results from the 15 day analysis of 
the current data. During this period the tidal modulation accounted 


for 21.35 per cent of the total fluctuations, with the remainder being 








Table 6. 


Statistical Results and Energy Calculations for 


the 15 Day Period. 
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TABLE 6 


Statistical results and energy 


calculations for the 15 day 


period 

mumpeemicman Gurremte Velocity . . . . .« « « » » « » » 167.05 cm/sec 
2 Standard Deviation ..... cm/sec 
oe Average Fluctuation... . per cent 
4. Total Current Variance. . Za 0 cm*/sec* 
5. Predicted Tidal Current Variance. cm*/sec* 
S: Fluctuations Due to Tidal Components. per cent 
ae Percent of Item 6 Due to Major Diurnal 

Tidal Components (K, and 0j) . 71 per cent 
oe Percent of Item 6 Due to Major Semi-Diurnal 

Tidal Components (My and S92) 20 per cent 
a Energy Contributed by Individual Tidal 


Constituents (cm¢/sec2) 
Ky ° ° ° ° ry ° ° ° e e e r e e e 


O71 


Q1 
Remainder Semi-Diurnal Constituents. 


Remainder Diurnal Constituents 


16. 


US) 


30 


40 
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fluctuations of an apparent non-periodic nature. 
B. COMPARISON OF RESULTS 

The analysis of the surface current data reveals large diurnal 
tidal modulation of the Florida Current. This is in agreement with 
the tidal analysis performed on the transport of the Florida Current 
by Richardson and Schmitz (op. cit.), and tidal analysis of transport 
fluctuations from studies of electric potential measurements by 
Wertheim (op cit.). It is further in agreement with Project MIMI's 
results where underwater acoustic signals transmitted across the 
Straits of Florida over long periods of time have shown prominent 
diurnal phase changes. (Steinberg and Birdsall, 1966), (Clark and 
Yarnall, 1967). 

Recently Zetler calculated the amplitude of the Kj tidal current 
in the Florida Straits required to conform to observations of the 
oscillating diurnal tidal water transport in the Gulf of Mexico 
(Zetler, 1968). The K, amplitude of 0.11 knots found from the har- 
monic analysis of the Monster Buoy data is in remarkable close agree- 
ment to his calculated value of 0.12 knots. 

The following is a comparison of the phase angles of the major 
diurnal constituents of the tidal current with the phase angles for 
the tide at the Patrick Air Force Base and Miami Beach shore stations. 
No harmonic constants are available for tide at the Monster Buoy 


Fate tiude. 
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TIDE LOCATION LATETUDE PHASE ( G) 
KL OL 

Patrick AFB 28°14' N 203 207 

Miami Beach 25) G6. oN 245 ZO7 


TIDAL CURRENT 


Monster Buoy 26-01 " aN 284 22 


If a wave is progressive, the difference in phase between tidal 
current and the resultant tide should be zero degrees at any location, 
whereas in a standing wave situation there should be a 90° difference. 
As can be seen, the differences found were inconsistent and probably 
indicate some combination of both a progressive and standing wave. 
Since, cotidal lines bunch up near a node and tidal constants at the 
latitude of the current observations are not available, the above 
comparison is of marginal suitability. No analysis of tidal currents 
at any of the shore stations is available for comparison. 

Zetler (1968) concluded, after consideration of the available 
known Kj and 0] phase angles at shore stations along either side of 
the Straits of Florida, that there is a strong indication of a longi- 
tudinal standing wave situation for the major diurnal tidal components 
in the Straits of Florida, with a node close to the latitude of Miami. 
The large amplitudes of the diurnal constituents of the tidal modu- 
lation of the Florida Current at the Hollywood latitude are in agree- 
ment with this concept as the tidal current related to this wave 


should be maximum at the node. 
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V. SUMMARY AND CONCLUSIONS 


During the latter portion of 1965 direct surface current measure- 
ments of the Florida Current were taken at hourly intervals over an 
approximate period of one month from the General Dynamics Monster 
Buoy anchored in the Straits of Florida at the Hollywood latitude. 

The data from this period presented the opportunity of conducting the 
first Florida Current tidal analysis from direct surface current 
measurements of sufficient duration to provide meaningful results. 

Harmonic analysis of data from 15 complete days of this period 
confirms that the resultant tidal influence on the Florida Current 
surface flow does not conform to the usual Atlantic coastal tidal 
configuration, being instead transitional between the semi-diurnal 
Atlantic tide and the diurnal Gulf of Mexico tide. The tidal coupling 
between the Gulf of Mexico and the Atlantic Ocean with pronounced 
diurnal features can only at present be explained by a heretofore 
overlooked longitudinal diurnal standing wave in the Florida Straits. 
The obtainment and subsequent analysis of tide observations at ad- 
ditional points along the lower third of the east coast of Florida 
should confirm the presence of this diurnal standing wave. 

A fluctuation of 10 per cent of the mean surface current is given 
as representative of the Florida Current. Slightly more than one- 
fifth of this modulation is attributed to tidal influence, being the 
apparent major periodic modulating force of the Florida Current. 


It is recommended that further work be continued in this field. 
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The obtainment of tidal current data in the Straits of Florida from 

a deep subsurface buoy would be particularly invaluable for analysis. 
Additionally, the placement of another data gathering platform in the 
Florida Current at about the Miami latitude for long term fluctuation 


study of the Florida Current would be well worth the expense. 
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APPENDIX A: 


DERIVATION OF THE ENERGY FORMULA OF 


A HARMONIC TIDAL COMPONENT 
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by: 


where 


Bye, 


If any harmonic component of the mean current speed is represented 


Vv 


W 


a,Sinwt 48 b ,Coswt 


= ba is the angular speed of the motion, then the 


energy of the periodic motion may be represented by: 


Pam + K.E. = 1 § vdt 
aaa 








T T 
=1 ag Sin-wtdt opie by “COS*wtdt se a, b,SinwtCOSwtdt 
it T u 
0 0 
T A u 2 ae 
= ar? a Sau abe 5+ ae + and, Sin re 
7 4w 0 TF 2 4w a 2w 0 
= a,2 + b,2 - a,?Sin2wT + b,Sin2wT + anb,Sin2wT 
= Zz — ae ar 
= 1 (a;+b,*) - (aq?-b,2) (Sin2wt) + (anbp) To 
2 a: (4wt ) 
=1 (a.2tb,.2) - (a.2-b,*) (1 b in@ 
LS kek fe 2a) 7 Ft (an n) Sin (ae) (x) 
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APPENDIX B: 


FLORIDA CURRENT DATA 


AND 


ENVIRONMENTAL DATA FOR 


THE 15 DAY PERIOD 
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TIME CURRENT SPEED WIND SPEED WIND 

DATE (LOCAL) (kts) (mph) DIRECTION 
11/21/65 0000 Bl 7 190° 
u 0113 3.4 6 185° 
i 0203 3.6 4 120° 
’ 0258 B..5 i ie O- 
0401 a. 5 3 120° 
‘ 0500 8.5 4. 100° 
m 0600 3.6 7 100° 
tk 0657 3.4 9 1G? 
i! 0800 3.4 10 100° 
m 0902 3.4 4 120° 
i 1000 oa4 4 GO? 
1100 Bt 5 5 Zo" 
: 1200 Ses 2 125° 
: 1300 Bed 2 120° 
y 1400 3.9 D 200° 
ut 1503 6.4 3 200° 
i 1600 4.0 4, 180° 
1700 oars 3 190° 
UL 1800 3.6 3 190° 
i 1900 3.0 d, 220° 
2001 3.2 9 200° 
A 2110 eee 8 200° 
uw 2206 Se 8 2R0- 
2 2303 3.3 10 220° 
11/22/65 0002 3.8 ia 260° 
" 0101 3.4 ll 250° 
u 0206 nD ik EQ? 
i 0301 Shae) ly, 240° 
0401 Beal 6 250° 
i 0502 2uell 16 260° 
0604 oD 14 270° 
. 0702 Bee 12 250" 
i 0800 ane 10 300° 
” 0920 cal 12 260° 
" 1000 Be 14 280° 
a 1100 3.2 10 290° 
i 1200 3.0 13 260° 
1301 Ba 13 260° 
5 1400 B54 is 272° 
0 1500 Be 14 2702 
iH 1600 ews 17 325° 
u. 1700 2.8 18 3007 
" 1800 3.5 19 Bao 
u 1900 Bae 421 3250 
‘ 2002 3.2 29 330° 
: 2104 353 16 3a0" 
: 2204 358 12 330° 
4 2303 ah 15 3256 
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TIME CURRENT SPEED WIND SPEED WIND 

DATE (LOCAL) (kts ) (mph ) DIRECTION 
ay 28/65 0000 363 2 330” 
: 0100 3.4 15 320" 
- 0200 3.2 19 320” 
= 0259 3.4 20 320° 
: 0402 3.4 18 320° 
: 0500 Diag 19 320° 
: 0603 ee 20 325° 
: 0700 3.2 20 340° 
a 0800 B20 i 340° 
- 0900 Lee, 12 020° 
e 1000 Zi 1y2 040° 
‘ 1100 Die} 12 030° 
. 1200 2.6 mE 060° 
‘ 1300 250 iz 060° 
- 1400 2.8 10 060° 
i‘ 1500 Dest 5 070° 
- 1600 a9 8 070° 
. 1701 2.8 8 080° 
- 1800 2.6 10 068° 
? 1900 a.2 16 IO” 
i 2001 a0 14 130° 
. 2100 3.0 17 130° 
2200 3.4 1 130° 
a 2302 San2 10 120° 
11/24/65 0002 2.2 Ae TO 
- 0110 31672 hs 20° 
fs 0220 a) 11 120 - 
7 0300 3.4 10 120° 
a 0400 vee 12 120" 
2 0500 BS: 15 120” 
- 0558 3.4 15 Og 
ii 0700 B15 5. INS 115° 
0800 Bees, 12 eo) 
" 0900 Bee 10 enial 
7 1000 ote 11 0: 
: 1100 om 12 be aye 
‘ 1203 a9) 15 12s 
‘ 1300 3.4 12 140° 
1400 oe ES he 2 
s 1500 Jeo 10 120° 
2 1600 3.4 ll 140° 
- 1700 oat LS 140° 
" 1800 eel 15 140° 
x 1900 oa 16 140° 
ui 2003 3.4 14 150° 
a 2104 325 10 [Sor 
‘ 2202 3.6 10 Ole 
2304 356 i 180° 
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TIME CURRENT SPEED WIND SPEED WIND 
DATE (LOCAL) (kts) (mph) DIRECTION 
11/25/65 0000 S65 20 170° 

is 0100 3.6 18 180° 
u 0201 5.8 18 180° 
: 0300 348 20 180° 
i 0400 3.6 19 200° 
is 0458 Sead 18 180° 
a 0600 3.8 14 200° 
u) 0701 3.8 15 200° 
uh 0800 Sa 15 188° 
Z 0900 3.6 10 140° 
i 1000 Sad 16 060° 
" 1100 &.5 20 195° 
i! 1200 3.4 1 060° 
: 1300 3.5 7 200° 
" 1400 oe. 4 060° 
am 1500 See 5 060° 
‘i 1600 2.2 5 Oy 5° 
i 1700 3.2 9 120° 
HL 1800 5.4 10 160° 
mn 1900 ce 9 140° 
i 2001 3.3 i 160° 
uy 2100 5.4 9 130° 
ts 2201 6.4 1? 150" 
As 2300 3.6 10 150° 
11/26/65 0000 305 2 kyo 
m 0103 3.4 10 iSo° 
i 0200 3.4 8 180° 
tt 0301 oa5 7 130° 
2 0401 Bn 5 10 15° 
“ 0501 3.6 6 160° 
u 0605 3.6 6 160° 
iy 0700 e058 9 140° 
1 1800 3.6 7 LoS 
. 0900 3.6 7 10° 
1000 3.4 10 lis 0° 
i 1100 323 9 Woh 
My 1200 3.0 10 140° 
i 1300 Bll 9 130° 
: 1400 $.0 14 140° 
1500 FES, 15 140° 
" 1600 3.0 We 170° 
n 1700 3.0 12 170? 
ie 1800 al 18 160° 
n 1857 3.0 19 We 
ui 2000 3.0 17 Lesh 
ss 2100 gy 18 180° 
2200 3.3 16 180° 
M 2257 BG 20 180° 
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TIME CURRENT SPEED WIND SPEED WIND 
DATE (LOCAL) (kts) ___ (mph) DIRECTION 
11/27/65 0000 3.5 20 LS Oba 

" 0058 3.6 18 190° 
" 0200 a65 18 190° 
tt 0300 3.6 19 220° 
" 0403 3.6 2 a0" 
" 0505 Suab 1 190° 
‘ 0615 3.6 2 Oe 
" 1700 3.6 5 iLO 
it 0810 3.7 9 Done 
J: 0900 Be 10 180° 
" 1000 3.5 10 160° 
i 1100 3.6 i I 70 
ge 1200 8.8 12 185° 
tt 1300 3.2 14 185° 
“ 1400 se | 10 200" 
ef 1500 3.8 6 220° 
4 1600 Bi 2 4, 200° 
p 1700 Se 7 200° 
st 1800 3.4 7 2U5° 
a 1900 Beh 8 200° 
rn 2000 8.8 8 200° 
i 2100 Ba 3 8 210° 
" 2200 3.2 8 210° 
i 2300 Bi 8 230° 
11/28/65 0000 3.8 9 250° 
‘i 0100 B43 9 270° 
u 0200 3.5 10 280° 
i 0300 3.6 10 330° 
M 0400 3.6 10 340° 
u 0500 3.8 Lil 340° 
A 0600 3.8 ll 340° 
i 0700 3.8 10 345° 
nt 0800 3.6 10 340° 
i 0900 3.8 10 340° 
1000 exc 10 350° 
ef 1100 oun 9 345° 
i 1200 Br? 8 330° 
i" 1300 3.0 5 340° 
ue 1400 3.1 3 268° 
5 1500 3.1 3 280° 
z 1600 3.1 4, 270° 
D 1700 Bel 6 275° 
et 1800 ae: 8 280° 
J. 1900 Be 10 260° 
A 2000 3.3 7 280° 
A 2100 oe 8 300° 
i 2200 Bek 12 220° 
: 2300 3.6 16 340° 





DATE 


11/29/65 


TIVE 
(LOCAL) 


0000 
0100 
0200 
0300 
0400 
0500 
0600 
0700 
0800 
0900 
1000 
1100 
1200 
1300 
1400 
1500 
1600 
1700 
1800 
1900 
2000 
2100 
2200 
2300 
0000 
0100 
0200 
0300 
0400 
0500 
0600 
0700 
0800 
0900 
1000 
1100 
1200 
1300 
1400 
1500 
1600 
1700 
1800 
1900 
2000 
2100 
2200 
2300 


WN WN NNW WW WWW WWW WWW WWW WWW WW WW WWW WW WW WWW WWW WWW WWW 
DOU NM DPOOOPEP REP RULWHIAKRUU A ~~ +1 GW MED OF em UO Rigo o Locate Cec. 


CURRENT SPEED 


(kts ) 


WIND SPEED 


(mph ) 


WIND 


DIRECTION 


040° 
CG: 
000° 
350" 
Sa: 
000° 
000° 
350 
010° 
O20. 
000° 
020° 
020° 
020° 
150° 
220° 
255 
240° 
ys 
280° 
280° 
280 
290) 
300° 
Sor 
320° 
320° 
000° 
030° 
Os oh 
060° 
030° 
035° 
350" 
000° 
005° 
000° 
350° 
050° 
050° 
060° 
050° 
060° 
060° 
060° 
060° 
060° 
050° 
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TIME CURRENT SPEED WIND SPEED WIND 
DATE (LOCAL) (kts) (mph) DIRECTION 
12/10/65 0000 S70 22 050° 
" 0100 2.8 25 060° 
q 0200 3.0 24 040° 
i 0300 2.8 25 058° 
ut 0400 2.8 25 058° 
p 0500 3.0 25 058° 
iL 0600 3.2 27 056° 
a 0700 3.2 27 055° 
i 0800 2.9 22 058° 
f 0900 304 25 060° 
e 1000 S38 22 060° 
M 1100 = al 20 060° 
7 1200 320 Ail 060° 
. 1300 2.7 20 060° 
L 1400 30 20 060° 
i 1500 go 19 060° 
Z 1600 3.0 20 080° 
1700 2.9 16 100° 
: 1800 3.0 20 eo? 
1900 2.8 20 090° 
" 2000 340 26 age 
, 2100 2.8 22 0° 
a 2200 370 22 Lier 
i 2300 3.0 20 120° 
2/2/65 0005 2.9 25 io 
n 0100 2.9 26 20° 
" 0200 2.9 22 140° 
4 0300 3.0 20 ler 
y 0400 2.9 20 LO 
et 0507 ee: 25 ie 
14 0602 728 22 25° 
: 0701 256 20 10° 
i 0800 27 25 iO? 
i‘ 0900 2x6 20 140° 
iM 1000 228 25 160° 
4 1100 2.8 19 ASO 
iu 1200 3.0 20 150° 
a 1300 2.6 19 160° 
n 1400 2.8 20 Ls 
" 1500 2.6 16 1502 
a 1600 2.6 19 160° 
y 1700 Dies 18 180° 
iN 1800 2.6 eS 180° 
H 1900 228 20 180° 
d, 2000 DS 18 160° 
4: 2100 su, 20 180° 
o! 2200 Cal 20 190° 
" 2300 3.0 22 160° 





DATE 


65 


TRAE 


(LOCAL) 


0000 
0100 
0200 
0300 
0400 
0500 
0600 
0800 
0800 
0900 
1000 
1100 
1200 
1300 
1400 
1500 
1600 
1700 
1800 
1900 
2000 
2100 
2200 
2300 
0000 
0100 
0200 
0300 
0400 
0500 
0600 
0700 
0800 
0900 
1000 
1100 
1200 
1300 
1400 
1500 
1600 
1700 
1800 
1900 
2000 
2100 
2200 
2300 


CURRENT SPEED 


(kts) 


Wo W WD WD WH WH WH WH WH W WD WH W WH WH WH Ww WWD W DW WW WWW WWWWWNWWNNMWDNHWNNNNN DN DH WH 


PERWWNWNHNPHENNOPFP FP OOOOVDOONWWNNODOHKPOODOUYNFWOOWMDNDBDIDRDDWO 


WIND SPEED 
(mph ) 


WIND 


DIRECTION 


160° 
150° 
140° 
70° 
190° 
180° 
150° 
io" 
160° 
160° 
170° 
170° 
GO - 
130° 
1ée 
180° 
Sia 
180° 
180° 
ae 
190° 
LOO” 
£90" 
Lesh 
200° 
200° 
200° 
200° 
200° 
200° 
200° 
200° 
180° 
200° 
190° 
220° 
240° 
2204 
230° 
240° 
240° 
280° 
OC0° 
065° 
060° 
070° 
060° 
060° 
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TIME CURRENT SPEED WIND SPEED WIND 
DATE (LOCAL) ___(kts) (mph) DIRECTION 

12/5/65 0000 3.5 18 080° 
0051 a4 15 100° 
i 0200 3.4 18 090° 
J 0300 one 15 090° 
i 0400 py 15 060° 
‘ 0500 ow 12 152 
u 0600 ep Le 100° 
u 0700 3.0 17 059° 
i 0800 eal 12 080° 
0900 Bel 16 058° 
i 1000 an iW 060° 
i" 1100 Sa5 ll 060° 
u 1200 3.6 8 LEO 
j 1300 3.6 12 Leer 
: 1400 3.5 ll LAO 
’ 1500 Bal 12 150° 
‘ 1600 3.6 10 2 5e 
s 1700 Sad 10 165° 
‘ 1800 3.3 6 200° 
p 1900 3.4 8 180° 
u 2000 3.6 5 ioe 
i 2100 ce: 6 Se 
2 2200 4.0 3 IO) 
uy 2300 4.0 5 MEY 








VITA 


Lt. John Alan Smith, USN, was born in Norwood, Massachusetts, 
on December 5, 1937. His parents were Arthur L. Smith and Elinor 
Smith. He received his elementary education in the Walpole public 
school system at Walpole, Massachusetts, and his secondary education 
at Fort Lauderdale High School, Fort Lauderdale, Florida. 

In September 1955, he enlisted in the United States Navy, and 
served in the USS NAUTILUS (SSN571) and attended the U.S. Naval 
Academy Preparatory School at Bainbridge, Maryland from September 
1956 to June 1957. In July, 1957, he entered the U.S. Naval Academy, 
Annapolis, Maryland. Upon graduation in June, 1961, with a B.S., he 
was commissioned an Ensign in the U.S. Navy. Subsequent Naval Service 
included five years in the Destroyer Force of the Atlantic Fleet with 
attendance at Fleet ASW School at Key West, Florida, and Fleet 
Destroyer School at Newport, Rhode Island. 

He was admitted to the Graduate School of the University of 
Miami in September, 1966. He was granted the degree of Master of 


Science in July, 1968 


Permanent address: 1911 North Sumac Drive 
Jaynesville, Wisconsin 
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